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Abstract 


Using superposition and numerical approximations of an analytical expression for the magnetic 
field generated by a finite solenoid, we show that the magnetic field external to parallel identical 
solenoids can be nearly uniform and substantial, even when the solenoids have lengths that are 
large compared to their radii. We study two arrangements of solenoids—a ring of parallel solenoids 
whose surfaces are tangent to a common cylindrical surface and to nearest neighbours, and a large 
finite hexagonal array of parallel solenoids—and summarize how the magnitude and uniformity of 
the resultant external field depend on the solenoid length and distances between solenoids. We 
also report some novel results about single solenoids, e.g., that the energy stored in the internal 
magnetic field exceeds the energy stored in the spatially infinite external magnetic field for even 
short solenoids. These results should be broadly interesting to undergraduates learning about 
electricity and magnetism as novel examples of superposition based on a familiar source of magnetic 
fields. 
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I. INTRODUCTION 


In introductory physics courses on electricity and magnetism 1 and in upper-level un¬ 
dergraduate courses on electrodynamics^, students learn a striking fact which is that the 
magnetic held B generated by a sufficiently long cylindrical solenoid is highly uniform inside 
and of small magnitude outside. An implication of this fact is that the external magnetic 
held generated by a solenoid can usually be ignored as small. 

But this well-known result raises some questions that are often not discussed in under¬ 
graduate courses. One question is how long must a solenoid actually be for the external 
magnetic held to be smaller than some specihed value, at different points in space? And 
might there be circumstances in which the external magnetic held of several solenoids could 
combine to be become substantial in magnitude, even for long solenoids? If so, might there 
also be configurations of solenoids such that their external magnetic held is also approxi¬ 
mately uniform, so that a single solenoid or a pair of Helmholtz coils 2 is not the only way 
to generate an approximately uniform magnetic held in some region of space? 



FIG. 1: Perspective drawings of two geometric arrangements of parallel identical solenoids, 
a) A ring of 5 solenoids, each of which is tangent to nearest neighbors and to a common 
inner cylindrical surface. The direction of the constant current density t on the surface of 
the solenoids is indicated by the arrows, b) A section of an extended hexagonal array of 

parallel solenoids. 

In this paper, we investigate these questions for two geometric arrangements of identical 
parallel solenoids (see Fig. |T]) , a ring of solenoids each of which is tangent to nearest neighbors 
and to a common inner cylindrical surface, and an extended hexagonal array. For both 
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cases, we find that there is a range of parameters (spacing of the solenoids and length of the 
solenoids for a given current density) for which the external magnetic field is substantial in 
magnitude, i.e., comparable to the magnitude of the field at the center of a long solenoid. 
Further, for a large hexagonal array of solenoids, the magnetic field between the solenoids 
is highly uniform and decreases slowly in magnitude with increasing solenoid length. 

As part of our analysis, we also discuss some calculations and insights regarding the 
magnetic field of a single finite solenoid. For example, we calculate how the external mag¬ 
netic field of a single solenoid decreases with increasing solenoid length for points ranging 
between close to and far from the solenoid’s surface and find that a simple analytical ex¬ 
pression derived for the far-held of a solenoid on its symmetry plane (see Eq. (J6]) below) 
gives a surprisingly accurate approximation of the magnetic held even close to the surface 
of the solenoid, at least for solenoids whose lengths exceed about four radii. This result 
unifies several previous discussions. We also show that the energy associated with the inter¬ 
nal magnetic held of a solenoid exceeds the energy associated with the infinitely extended 
external magnetic held once the solenoid length exceeds about the radius of the solenoid, so 
comparing the internal and external magnetic energies does not lead to a new length scale 
that could be used to distinguish long from short solenoids (see Fig. [5] below). 

While most of our results are new, parts of our Section [TJ overlap with some previous 
papers. Brown and Flax- 4 discuss a way to calculate the magnetic held of a thick solenoid 
starting from the same integrals we use- 11 ’ for a zero-thickness solenoid, and Derby and 
Olbert- has a closely related discussion that discusses a simple code to approximate these 
integrals using elliptic functions. Farley and Priced show that the magnetic held strength 
just outside and on the midplane of a long finite solenoid of arbitrary constant cross-section 
falls oh as L~ 2 with increasing L but these authors do not discuss, as we do for solenoids 
of circular cross section, how the magnetic held behaves for short solenoids or for points 
that are at intermediate distances from the solenoid’s surface. Muniz et al.- used Taylor 
expansions of the magnetic scalar potential to calculate the off-axis internal magnetic held 
of a finite solenoid, which complements our direct approximation of two one-dimensional 
integrals using numerical integration, although our method works for any point in space. 
However, these earlier papers do not mention some of our single solenoid results such as how 
the external magnetic energy compares with the internal magnetic energy, or how there is 
a simple expression that provides an accurate approximation of the magnetic held on the 
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midplane of a solenoid, for all points, not just close to the solenoid’s surface. We also do not 
know of previous discussions about the magnetic field created by multiple parallel solenoids, 
which is the key contribution of this paper. 

Our results should be of broad interest to undergraduate physics students and to in¬ 
structors of undergraduate physics courses as a moderately more complicated example of 
superposition of magnetic fields, based on an example—a finite solenoid—that students have 
already learned about. The example we discuss below for a ring of parallel solenoids shows 
that just having a high symmetry of sources (many identical solenoids tangent to an inner 
cylinder) is not enough to guarantee uniformity of the external magnetic field. The second 
example of a large hexagonal array of parallel solenoids also shows that the external magnetic 
field of solenoids can be substantial, even for solenoids whose lengths are large compared to 
their radii. Finally, this paper provides a useful example to share with students of how soft¬ 
ware like Mathematical, Maple^, and Matlab^ makes it easy for undergraduates to study 
numerically and visually superpositions of magnetic sources that are complicated by their 
number or geometry. 

The rest of this paper is organized as follows. In the next section, we discuss a previously 
published analytical result^- for the magnetic field generated by an idealized cylindrical 
continuous solenoid and some insights about this field based on numerical studies of the 
analytical solution. In Section UTT1 we use the results of Section [Til to discuss the properties 
of the magnetic field generated by a ring of identical parallel continuous solenoids, while in 
Section m we discuss the properties of the magnetic field generated by a large hexagonal 
grid of identical parallel continuous solenoids. After summarizing key points in Section [V] we 
discuss in Appendix lAl how various numerical calculations were validated, and in Appendix iBl 
how the radial components of the magnetic fields from different solenoids were combined. 


II. THE MAGNETIC FIELD OF A SINGLE FINITE CONTINUOUS SOLENOID 

The starting point for our study of the external magnetic field created by superimposing 
the fields of several parallel finite solenoids is an analytical expressions^ detailed below, 
for the magnetic field B generated by a continuous solenoid, by which which we mean a 
cylindrical surface of radius a and length L such that a spatially uniform time-independent 
one-dimensional current density i (with units of amperes per meter) flows azimuthally around 
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the surface. 

In this section, we discuss some physical properties of this analytical expression, for 
example we characterize how the external magnetic field depends on its length L, and 
confirm and extend an analytical result below that the external magnetic field at some point 
in space first increases and then decreases with increasing L , with an asymptotic function 
behavior of L~ 2 . We also compare the energy associated with the external magnetic field 
with the energy associated with the internal magnetic held as a function of L/a and find 
that the two energies become comparable when L « a, so that there is not an interesting 
new length scale for solenoids related to the magnetic energy. 

A continuous solenoid is an idealization of real solenoids that are helically wound with 
wires of some finite thickness. The continuous solenoid has the advantage over real solenoids 
of requiring just three parameters to specify—a radius a, length L, and uniform current 
density t —and has the further advantage that the analytical expression for its magnetic 
held is easy to evaluate numerically, accurately and quickly, at any point in space. The 
continuous solenoid is somewhat unphysical in that the zero-thickness of the surface causes 
the radial component of the magnetic held to diverge in magnitude at the top and bottom 
edges of the solenoid’s surface (see Fig. [3](a)), while no such divergence occurs for physical 
solenoids of finite thickness. However, we show below in this section that this divergence is 
not important since it still leads to a finite amount of energy stored in the magnetic held, 
and that the region where B r has a large magnitude has a small spatial extent (compared 
to the solenoid’s radius) so that the magnetic held of the continuous solenoid provides 
quantitatively useful insights about the magnetic held of physical solenoids. 

To describe the axisynnnetric magnetic held of a continuous solenoid of radius a, length L , 
and current density t, we introduce a cylindrical coordinate system (r, 6, z ) such that r = 0 
defines the solenoid’s axis and the plane z — 0 bisects the solenoid (so the ends of the solenoid 
lie at the coordinates z = ±L/2). For this coordinate system, the current density t hows 
counter-clockwise, so that the solenoid’s internal magnetic held B points in the positive z 
direction. From the azimuthal symmetry of the problem, we deduce that Bg = 0, and 
furthermore that B z and B r do not depend on 6, so that the magnetic held has the form 

B = B r (r, z)r + B z (r, z)z. (1) 

Calculations then show^— that the magnetic held at a point (r, z) in space has components 
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given in terms of the following one-dimensional integrals: 
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The lengths are dehned by 

^± = z±-, (4) 

and the notation [/(£)]^ + is an abbreviated way of writing /(£+) — /(£_). 

Appendix [A] gives details of how we approximated Eqs. ([2]) and © numerically for given 
values of r, z, a, L , and t using Mathematical. This appendix also summarizes how we 
validated the numerical approximations of these integrals, and results obtained for super¬ 
positions of solenoids. We found that Mathematica’s numerical approximations to Eqs. (}2]) 
and (|2]) had a relative accuracy that exceeded eight significant digits, except for the compo¬ 
nent B r near the coordinates (r,z) = (a, ±L/2), where this component diverges to infinity. 
This accuracy was more than adequate for our calculations. 


Z 



FIG. 2: Vector plot of the axisymmetric magnetic held B(r, z) generated by a finite 
continuous solenoid, Eqs. (DO)-©, for radius a — 1, length L — 2, and current density i — 1. 
Vectors are shown only for r > 0 so only the lower half of the solenoid (gray region) and 
half of the magnetic held are shown. The internal held is approximately uniform despite 
the solenoid’s short length. Even for this short length, the interior magnetic held has 
attained a magnitude B ~ 0.8-E> 0 at the solenoid’s center r = z = 0, where B 0 is dehned in 

Eq. ©. 
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Fig. [2] shows a vector-field plot of the magnetic field B generated by Eqs. ([2]) and (j3J) for 
a solenoid of radius a — 1, length L — 2, and current density i — 1. Even for this short 
length, the interior magnetic field is approximately uniform and the magnetic field at the 
solenoid’s center r = z = 0 has attained a magnitude B rs 0.8-B 0 that nearly equals the 
magnitude 

-Bo — hob (5) 

of the uniform internal magnetic field of an infinitely long solenoid with the same current 
density i. We note that Eq. (J5]) is the familiar formula^* B = Honl for the magnetic field 
magnitude inside a long solenoid that is uniformly wound with n wires per unit length and 
with each wire carrying a current I. The product nl has units of current per length and 
corresponds to the current density i in the limit n — >■ oo, / — y 0 with nl finite. 


(a) (b) 



FIG. 3: The magnetic field B(r, z) of a finite continuous solenoid has two discontinuities, 
here shown for parameters a = 1, length L — 4, and unit current density t = 1. (a) The 
radial component B r (r, z ) diverges for z = ±L/2 at r = a. Here we show four curves 
of B z (r,z) for z = 0.8(L/2), 0.95(L/2), 0.99(L/2), and L/ 2 (from lowest to highest curves). 
The curve for z = L/2 (and similar curve for z = —L/2) shows a divergence of B r at r = a. 
One can show from Eq. (2) that this is a logarithmic divergence of the form ln(|r — a|) for 
r sufficiently close to a and for 0 = ±L/2. (b) The axial component B z (r, z) has a jump 
discontinuity for —L/2 < z < L/2 at r = a that represents the change in direction and 
magnitude of the magnetic field as one progresses radially from just inside to just outside 
the surface of the solenoid on a line of constant z. 

It is difficult to see from Fig. [2] that the magnetic field of a finite continuous solenoid has 
two kinds of discontinuities. Fig. |3](a) shows how the radial component B r (r, z) varies with r 
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for fixed z, for several different values of z that approach a singularity at r = a that occurs 
for z = ±L/2. (So B r diverges everywhere on the two ends of the solenoid at z — LL/2, 
which is again an artifact of our approximating the solenoid as a zero-thickness surface; this 
singularity does not occur for a solenoid of finite thickness with a corresponding volumetric 
current density J .) By setting z — L/ 2 in Eq. (2), one can evaluate the integral analytically 
over the range 0 < r < a — e and then find that the integral has a term proportional to ln(e) 
and so diverges logarithmically in the limit e —>■ 0. This will imply that the total energy 
associated with the magnetic held (Eq. (7) below) is a finite quantity. 

The surface plot of the axial component B z (r,z) in Fig. El(b) shows a more familiar 
discontinuity, that the ^-component B z of the magnetic held abruptly changes direction and 
magnitude as one passes radially from just inside to just outside the surface of the solenoid. 
(This is also apparent from Fig. [2] for r near r = 1 and z — 0, you can see how the magnetic 
held B changes direction and magnitude near the surface of the solenoid.) 

We next discuss how the external magnetic held of the finite continuous solenoid decreases 
with increasing length L for hxed radius a and for hxed current density i, an important issue 
for our discussion in later sections regarding how large can the magnitude of the external 
magnetic held be for multiple parallel continuous solenoids of some given length L. Here 
there is some prior analytical insight via Exercise 5.61 on page 254 of Griffith’s textbook on 
electro dynamics^, which states that the magnetic held on the solenoid’s midplane (z = 0) 
far from the axis (r ^ a) asymptotically has the form B = —B z, where the magnitude B is 

La 2 


B ps /iqL x 


4(r 2 + (L/2) 2 ) 


3/2 ‘ 


( 6 ) 


This is the leading-order term in an expansion in powers of the small quantity a/r and so 
its validity when a/r is not small, say close to the surface of the solenoid, when r ~ a, is 
not known. 

Equation ([6]) makes two predictions. First, for sufficiently large solenoid lengths L r, 
B~ oc L~ 2 so B z decays rather slowly (algebraically rather than exponentially) to zero with 
increasing L. We note that this 1/L 2 scaling is also how the magnetic held at the center of 
the solenoid, r = z = 0, converges to its infinite-length value Eq. (J5J) i.e. (B z — B 0 )/B 0 oc L -2 
for sufficiently large L. (This follows from the analytical expression for the magnetic held 
on the axis of a finite solenoid, see Eq. (1A8I) in Appendix [Aj ) Second, if this expression has 
validity for general values of r, Eq. ([6]) predicts that at a hxed radial distance r from the 
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axis, B does not decrease monotonically with increasing L, but first increases with a local 
maximum at L max = \/2r, and then asymptotically decays to zero as 1/L 2 . 

In Fig. U we compare these predictions with the exact analytical held of Eqs (J2]) and (J3]) 
for radial distances that are close to and far from the solenoid’s surface. Panel (a) shows that, 
for the representative fixed radial values r = 1.1a, 2a and 4a on the z = 0 symmetry plane, 
the ^-component of the magnetic held B z (r,L) indeed hrst increases and then decreases 
with increasing L as suggested by Eq. (]6jh even for relatively short solenoids for which 
L/a ^ 2. Surprisingly, the L-dependence of the analytical held at some hxed r is accurately 
described by Eq. ([6]) even for distances close to the solenoid (r > 2a). Panel (b) shows a 
complementary result, that even for modest hxed solenoid lengths (L > 4a), Eq. ([6]) also 
describes the radial dependence B z (r, 0) well, even for radii close to the outer surface of the 
solenoid. We conclude that Eq. (J6|) can be used as a quick and easy way to determine the 
magnitude of the magnetic held on the symmetry plane z = 0 of configurations of parallel 
solenoids without having to evaluate the integral Eq. ((3]). 

The analytical result Eq. (1A8h for B z on the axis of a finite continuous solenoid shows 
that solenoids whose lengths L > 4a are already “sufficiently long” in the sense that the 
magnetic held magnitude \B z [r = 0, z = 0)| at the solenoid’s center already exceeds 90% of 
its infinite-length value /i 0 l. It occurred to the authors that an alternative way to identify 
when a solenoid is “sufficiently long” would be to ask for what length L* (for a given radius a 
and given current density l) does the energy U ex t associated with the magnetic held 

U = f — — d 3 r = —— [ [ 2ttt ( B 2 + B 2 ) dr dz (7) 

J 2/i 0 2/i 0 J J 

external to the solenoid (the infinite region defined by r > a or \z\ > L/2), become less than 
half of the total energy U to t associated with the entire magnetic held (Eq. (J7J) evaluated over 
all of space, r > 0). By approximating the integral Eq. (J7|) numerically (see Appendix [All . 
we calculated and show in Fig. Othe ratio U ext /U to t, as a function of solenoid length L for 
parameters a = 1 and t = 1. This ratio falls below the value 1/2 for L « 0.69a so that, even 
for quite short solenoids, the energy associated with the magnetic held inside the solenoid 
exceeds the energy associated with the external held, even though the latter has infinite 
spatial extent. 
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FIG. 4: (a) Comparison of Griffiths’ result, Eq. (|6]h with the external magnetic held 
B z (r, z = 0) of a single finite continuous solenoid, normalized to the infinite solenoid axial 
value B 0 (Eq. (]5j)) versus the solenoid length L/a. Results are shown for three 
representative radial values, with the dotted lines showing Griffiths’ corresponding result. 

Note how B z first increases and then decreases with increasing L. For r = 4a, the 
analytical result is the same as Eq. ((6]) to the accuracy of the plot, (b) The asymptotic 
result Eq. (|6]) for LL(r, 0) derived for r a (dotted lines) accurately describes the radial 
dependence of the analytical held B z (r ) even for r close to the surface of the solenoid, 
provided that L > 4a The analytical result and numerical results are indistinguishable at 

the level of these plots for L/a > 7. 


III. THE MAGNETIC FIELD OF A RING OF IDENTICAL PARALLEL FINITE 
CONTINUOUS SOLENOIDS 

We are interested in the question of whether the magnetic held external to multiple 
parallel finite solenoids can ever be substantial in magnitude (comparable to the value B 0 = 
/io l on the axis of an infinite solenoid) and approximately uniform. In this section, we 
consider an arrangement of n > 3 parallel identical finite continuous solenoids, each with 
current density t = 1 and of length L, that are arranged in a ring as shown in Figs. [6] 
and IHa). Each solenoid is tangent to an inner common cylindrical surface of radius R and 
tangent to its neighbors. Because a substantial portion (nearly 50% in the case of many 
solenoids) of the external magnetic held of each solenoid passes through the common central 
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FIG. 5: Ratio of the energy U ext , Eq. 0. associated with the magnetic field external to 
the solenoid (r > a or \z\ > L/2) to the total magnetic energy U to t over all of space for a 
single continuous solenoid of fixed radius a = 1 and fixed current density t, as a function of 
solenoid length L. The energy associated with the magnetic field inside the solenoid 
already exceeds the energy associated with the external magnetic field for short solenoids 

such that L > 0.69a. 


cylindrical region, this ring geometry is a plausible candidate for producing an external 
magnetic field of substantial magnitude. We show in this section that this ring geometry 
can indeed produce a substantial magnetic held for solenoids that are long compared to their 
radii, but that the magnetic held in the common cylindrical region is generally not uniform. 
In the following section, Section IIV1 we discuss similar questions for a large finite hexagonal 
array of parallel finite continuous solenoids. 

To describe the external magnetic held in the central cylindrical region, we introduce 
a second cylindrical coordinate system (p, 0, z) such that p — 0 corresponds to the axis 
of the inner cylindrical surface, p = R corresponds to the cylindrical surface that all the 
solenoids are tangent to, and z = 0 is the common bisecting plane of the parallel solenoids. 
For n parallel solenoids that are all tangent to the inner cylindrical surface of radius R and 
tangent to neighbors, the tangency conditions imply that the common solenoid radius a n is 
given by 




s in (n/n) 

1 — sin(7r/n) 


n > 3, 


( 8 ) 


so a 3 ~ 6.5i?, a 4 ~ 2.4/2, a 5 ~ 1.4/2, a 6 = R, and a ~ ttR/ n becomes small for large n. 
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FIG. 6: Cross-section of a ring of n — 5 parallel continuous solenoids of length L — 2 and 
radius r « 1.4 that are tangent to neighboring solenoids and to a common inner cylindrical 

surface of radius R = 1, shown as a dotted ring with center C. The line segment CB 
denotes the location of the rectangular area, transverse to the plane of the figure, in which 
the magnetic held B is plotted in Fig. [7J For a perspective drawing, see Fig. Ufa). 

At any point (p, 0, z) inside the central cylindrical region p < R, the ^-component of 
the magnetic held B z (p,cj),z ) is obtained by directly adding the values B z (ri,z) from each 
solenoid in the ring at that point, where r t is the distance of the point ( p,(j>,z ) to the axis 
of the ith solenoid. Calculating the radial component B p (p, 0, z) is more involved since the 
radial direction r* with respect to the ith solenoid at (p, 0, z) varies with i. The details are 
given in Appendix iBl 

We now discuss our calculations of the magnetic held in the central cylindrical region 
for rings with varying number n of continuous solenoids and varying length L, for a cen¬ 
tral cylindrical region of hxed radius R — 1 and for all solenoids having the same current 
density i = 1. 

For a ring of n = 5 solenoids of radius a ~ 1.4 (see Eq. (jHJ)) and length L = 4a, Fig. [71(a) 
shows a vector plot of the magnetic held B(p, 0, z) inside the central cylindrical region 
0 < p < R along the plane 0 = 0, which corresponds to a rectangle that is perpendicular 
to and passes through the line segment CB in Fig [6](b) . We see that the magnetic held 
is approximately uniform in direction and magnitude with a magnitude B / Bq k, 0.3 that 
is about one third the magnitude of the magnetic held B$ = /i 0 t found on the axis of an 
infinitely long single solenoid with the same current density i. The magnetic held deviates 
substantially from that of a single solenoid (see Fig. [2]) in the planes z = ±L/2 rs ±2.8, 
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(a) (b) 

FIG. 7: External magnetic field produced by 5 parallel continuous solenoids tangent to an 
inner cylinder of radius R = 1, for solenoid length L = 4a ~ 5.6 where the solenoid 
radius a ~ 1.4 . (a) Vector plot of B (p, z) in the plane 0 = 0 (line CB of Fig. |6](b)) shows 
the direction and magnitude of the magnetic held inside the central cylindrical region. The 
lengths of the vector heads are proportional to the local magnitude B. (b) Surface plot of 
held magnitude B / Bq normalized to the magnetic held magnitude Bq inside an infinite 

single solenoid of the same current density t. 

which are indicated by the horizontal black lines in Fig. [71(a); the magnetic held is purely 
horizontal and of greatly decreased magnitude. 

Further insight about the structure of the central magnetic held is provided by panel (b) 
of Fig. [7] which shows a surface plot of the held magnitude B(p,z) for the same rectangle 
as panel (a). Panel (b) confirms that the magnetic held is approximately uniform inside the 
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FIG. 8: Plot of B z (p,<j),z = 0) inside the inner cylindrical surface of an arrangement of five 
solenoids each with L/a = 4 and t — 1, for p = 0.07, 0.29, 0.64 and 0.93. The magnetic held 
has high azimuthal symmetry, despite there being five discrete sources of magnetic held. 
There is only evidence of a hve-fold symmetry close to the edges of the solenoids. 


central cylindrical region (the magnitude B is approximately constant for \z\ < L/ 2) but 
becomes non-uniform near the boundary of one of the solenoids because of the divergence 
of B r near the ends of the solenoid (see Fig. 0(a)) . 

Fig. [8] shows that the magnetic held inside the cylindrical region is highly axisymmetric 
for z — 0. Despite the fact that the held is being generated by hve discrete solenoids, there 
is graphical evidence of a hve-fold symmetry only for r close to R amd for z — ±L/2, as 
shown in Fig [71(b). 

We found that the external magnetic held in the central cylindrical region becomes more 
nonuniform and decays radially in magnitude for larger numbers of solenoids (n > 4) so 
Fig. [7] is close to the best case in terms of achieving a uniform magnetic held in a ring 
geometry. As the number n of solenoids increases for hxed R and hxed L , the solenoid 
radii a n Eq. ([S]) decrease in size, the solenoids effectively become longer compared to their 
radius, and so the magnetic held external to each solenoid decreases in magnitude according 
to Eq. (J6J), being largest near the surfaces of the solenoids (p ~ 1) and decreasing towards 
the center of the common cylindrical region where p = 0. This is illustrated in Fig. [9] which 
shows how B z /B 0 (p, </> = 0, z = 0) varies on the bisecting plane z — 0 for increasing values 
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FIG. 9: Axial component B z (p , z = 0) on the bisection plane z — 0 for n = 3, 10, and 50 
solenoids, for fixed solenoid length L/a = 4 and fixed current density i. The increase in the 
number n of sources cannot compensate for the 1/L 2 decrease in magnetic held strength so 
the magnetic held in the central cylindrical region becomes small for large n. 


of solenoid number n, for a hxed length L/a = 4. 


IV. THE MAGNETIC FIELD OF AN HEXAGONAL ARRAY OF IDENTICAL 
PARALLEL CONTINUOUS SOLENOIDS 

The second conhguration of parallel finite solenoids that we explore is a large finite 
hexagonal lattice. Since we have shown in Fig [4] of Sec. [HI that the external magnetic held of 
a single solenoid can be as large as 0.25-B 0 on axis, just outside the solenoid, this geometry 
is a potential candidate for producing an external magnetic held of substantial magnitude 
(at least for specihc solenoid lengths and spacings). 

We define an hexagonal Bravais lattice 12 - such that the locations v of each solenoid axis 
in the z = 0 plane are given by the two-dimensional vectors 

v = s (mvi + nv 2 ) + v 3 , (9) 

where the coefficients m and n go over all possible integers, and where the positive parame¬ 
ter s denotes the lattice spacing (distance between two nearest neighbor solenoid axes). The 
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FIG. 10: (a) Horizontal cross-section through a portion of an hexagonal lattice of parallel 
continuous solenoids for lattice spacing s = 1 as defined by Eq [9] and [101 The common 
radius of the solenoids was chosen to be a = 1/4. Also shown are the directions of the 
coordinates x and y.(b) Perspective drawing of a section of the lattice. The black dashed 
rectangle corresponds to the region in which we compute the magnetic field shown in 

Figure [T2l 



basis vectors v x and v 2 are given by 

v x = i, and 



( 10 ) 


The vector v 3 , which here is (1/2) i + (l/2\/3)j, shifts the lattice relative to the origin of 
the coordinate system and places the origin at the centroid of the triangle formed by three 
nearest-neighbour solenoids. For a given lattice spacing s, the solenoid radius a must satisfy 
a < s/2 so that that the solenoid surfaces do not intersect. 

The array of solenoids is characterized by the lattice spacing s and the length-to-radius 
ratio L/a of each of the solenoids. We introduce an xyz-Cartesian coordinate system, where 
the origin ( x,y,z ) = (0,0,0) corresponds to the centroid of the triangle formed by the 
centers of three nearest-neighbour solenoids in a cell of the hexagonal lattice. Again, z — 0 
corresponds to the bisecting plane of the solenoids. 

We approximate the total magnetic field near the origin for an infinite period lattice was 
found by summing the fields from many solenoids in a large, finite, approximately circular 
region centered on (0,0,0). Fig. fill shows the component B z at (0,0,0) as a function of 
the radius of the area from which fields of individual solenoids were summed. Summing the 
fields from an area of the hexagonal array with radius 100s is sufficient to approximate the 
magnetic field due to the infinite array with a relative accuracy better than 0.05. 
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FIG. 11: Plot of the magnetic field component B z /Bq at the centroid of three solenoids in 
an hexagonal lattice as a function of the radius R of a circle centered at the centroid. For 
the largest radius R/s = 500 (not shown), there are ~ 750,000 solenoids enclosed within 
the circle. Note how over the range 30 < R < 500, B z varies only be about 10%, and 
furthermore appears to asymptote to a finite value. A circle of radius = 100s already 
approximates the infinite lattice case to better than 5%. 


The magnetic held at a given point ( x , y , z) was computed by individually summing 
the components of the magnetic held from the solenoids in the lattice from a region with 
radius 100s (see Appendix[B|). The held within each unit cell of the hexagonal lattice is again 
axisymmetric. Further details of the structure of the magnetic held within a single cell of 
the lattice are shown in Fig. [T2], which is a vector held plot of the magnetic held in the plane 
defined by the line section marked BC in Fig.^UJand with axial coordinates between z = ±L. 
The external magnetic held away from the edges of the solenoid is comparable in uniformity 
to the magnetic held inside the solenoid, with a magnitude B z = 0.6 Bq. This external 
magnetic held is quite uniform for most of z < \L/2\ , unlike the central magnetic held 
described in Sec. m 

We also found that the magnetic held B z (0, 0,0) decreases approximately algebraically 
as a function of the lattice size s for hxed a and L. Fig [13] shows how the magnetic 
held Fh(0,0,0) decreases in magnitude with increasing lattice size s for parameter values 
l — 1, a — 1/4, and L = 0.1,5,50, and 100. The decay of B z to 0 follows an approximate 
power law with a numerical exponent that approaches —2 with increasing solenoid length. 
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FIG. 12: Plot of the magnetic vector field B generated between solenoids in a hexagonal 
lattice in the black dashed rectangle shown in Fig. flOf bb with s = l,a = l/4, r = l, 
and L/a = 8, summed over solenoids in a circular area with radius R = 100. The location 
of the solenoid is indicated by a gray rectangle. Only the in-plane components of the 
magnetic field are non-zero. The magnetic field is highly uniform in the region external to 

solenoids, for \z\ < L/ 2. 


Furthermore, we find that the magnitude of the magnetic field for any given lattice spac¬ 
ing depends nonmonotonically on the length of the solenoid. However, this dependence is 
small: for fixed s, we observe in Fig. [13] that changing the lattice spacing by a factor of 20 
decreases B z /B 0 only by a factor of ~ 0.5. 

This slow decay of B z /B 0 towards zero as a function of L, for fixed s, can be seen more 
clearly in Fig. [HI which shows H z (0, 0, 0)/B 0 as a function of L/a for a = 1/4, s = 1, and 
l — 1. As expected from our discussion of the external magnetic field of the single solenoid, 
the decay of B z is nonmonotonic. However, after reaching a local maximum, B z decays 
linearly, with a numerically derived functional form B z = —0.0045L + 0.63, from which we 
extrapolate that B z decays to 0 only for a solenoid length of ~ 140. This linear decay of the 
external magnetic field of the single solenoid provides another example of a feature of the 


18 









FIG. 13: Plot of the ratio B z ( 0, 0, 0)/Bq as a function of the lattice spacing s for an 
hexagonal array of parallel identical solenoids of length L = 0.1, 5, 50 and 100, 
radius a = 1/4, and current density i — l. The decay of B z towards 0 follows a power law 
with a numerically derived exponent that approaches —2 for increasing L, consistent with 

Eq. ©. 

magnetic field that would have been difficult to predict without carrying out a numerical 
study. 


V. CONCLUSIONS 

In this paper, we have investigated using elementary methods that sophomore physics 
and engineering students can understand (numerical approximations to one-dimensional 
integrals and superposition by direct summation) whether the magnetic held external to 
some arrangement of identical parallel finite solenoids of radius a, length L , and uniform 
current density l can ever be substantial in magnitude and possibly uniform. Our motivation 
was to provide a deeper insight to the properties of the magnetic held external to a finite 
solenoid, which in many undergraduate physics textbook is ignored because its magnitude is 
known to become negligibly small for solenoids that are sufficiently long compared to their 
radii (although what is meant by “sufficiently long” is rarely discussed). 

We were also interested in the question of whether the supposedly small external magnetic 
held could ever be made substantial (compared to the magnetic held magnitude Eq. ([5]) 
inside an infinitely long solenoid of identical current density) and possibly highly uniform 
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FIG. 14: Plot of the ratio B z { 0, 0, $)/Bq as a function of the solenoid length L for an 
hexagonal array of solenoids. For L/a = 10, B z /Bq reaches a value of 0.58, as marked by 
an arrow. The magnetic field attains a local maximum at L — 2 of about half the field of 
a single finite solenoid, after which it decays slowly towards the infinite length limit 

of B z = 0. 

by combining the magnetic fields of many parallel identical finite solenoids. This question 
is not easy to answer, say by applying Ampere’s law to some small rectangular loop as is 
commonly done in undergraduate physics textbooks to a single infinite solenoid-^-, since one 
does not know in advance the symmetry of the magnetic held, e.g., whether it is everywhere 
parallel to the axes of the solenoids, and uniform in magnitude. 

We studied two configurations of identical parallel solenoids with high symmetry since 
high symmetry should favor uniformity of the external held. One configuration consisted of 
a ring of n parallel solenoids that each are tangent to nearest neighbors and to a common 
inner cylindrical surface (see Fig. [6]). The second configuration was a large, finite, roughly 
circular hexagonal array of parallel solenoids (see Fig. [10]), that we used to approximate the 
magnetic held generated by an infinite hexagonal array of parallel solenoids. For both cases, 
we investigated how the external magnetic held depended on the radius, length, and spacing 
of the solenoids, while the current density was kept hxed. 

For the ring geometry and hexagonal array, we found that the external magnetic held 
conld have a substantial magnitude, even for solenoids that are long compared to their radii. 
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But only for the magnetic field generated by a large hexagonal array did we find that the 
external magnetic field could also be highly uniform, at least in regions of space that lies 
between the two planes that contain the ends of the solenoids. 

The starting point for our analysis was an exact analytical expression, Eqs. (H])-Eq. (J3]) , 
in terms of two one-dimensional integrals for the axisymmetric magnetic field B produced 
by a single finite continuous solenoid. It turns out that these one-dimensional integrals can 
be rapidly and accurately approximated at all points in space (except for points near the 
edges r = a and z = ±L/2 of the solenoid, where the radial component B r diverges logarith¬ 
mically) to better than eight significant digits using simple and straightforward invocations 
of functions provided by the mathematics program Mathematical (see Fig. [15] as an exam¬ 
ple). The same Mathematica program provided easily used tools to add up the magnetic 
fields of many solenoids by superposition, and then visualize the resulting magnetic field. 

Thus this paper, in addition to providing new insights about the magnetic field generated 
by a single solenoid and by groups of solenoids, should provide a useful example to share 
with undergraduate physics students about how mathematics software environments like 
Mathematica can be used to explore physics problems that are substantially more challenging 
and more interesting than what can be analyzed by analytical methods. Based on our 
experience, we feel that such software should be more widely and systematically incorporated 
into undergraduate physics courses so that students learn how to use analytical, numerical, 
and visual ways to understand or explore a variety of physics problems. 


Appendix A: Validation of the numerical calculations 

The Mathematica code needed to carry out the calculations discussed in previous sections 
is modest, about four pages in length. But it took some effort to determine whether the 
results were correct and their accuracy. We summarize in this section some of the steps we 
took to validate the results. 

We used the Mathematica function NIntegrate to approximate the key one-dimensional 
integrals, Eqs. ([2D and ([2D, that give the cylindrical components B r and B z of the mag¬ 
netic field generated by a finite solenoid. Details about this function are available through 
the freely accessible URL reference.wolfram.com/language/ref/NIntegrate.htral, so 
we mention briefly that this function takes as its arguments some integrand in symbolic 
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form and the bounds of the definite integral (see Fig. [15]) plus many possible optional pa¬ 
rameters. (We used the default parameters and obtained acceptable accuracy as mentioned 
further below.) This function then uses an unspecified adaptive integration algorithm to 
numerically approximate the definite integral to some specified accuracy. 


bzSolenoid [ r_, z_, a_, L_ ] := Module[ 

{^p, ?n, 1 = 1., M° = 4. * n * 10' 7 } , 

= z + L / 2 ; 
ijn = z - L / 2 ; 

a J10 i 

- * NIntegrate 

2 n L 

(a-r Cos [a]) ^p 

(r 2 + a 2 - 2 a r Cos [a] ) (Jjp 2 + r 2 + a 2 - 2 a r Cos [a]) 0,5 (^n 2 + r 2 + a 2 - 2 a r Cos [a] ) 0-5 

, {a, 0, n} 


FIG. 15: Mathematica code that uses the Mathematica adaptive numerical integration 
function NIntegrate to approximate the magnetic field component B z {r, z,a, L) given by 
Eq. (|3]) at a point (r, z) in space, for a solenoid of radius a, length L, and current 

density i — 1. 


We tested the accuracy of the numerical values returned by NIntegrate for Eqs. ([2]) 
and ([3]) by comparing their values with analytically equivalent expressions that use entirely 
different numerical algorithms for their approximation. Researchers have shownr^ that 
Eqs. (J2]) and (J3]) can be expressed in terms of elliptic integrals^ as follows 
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where the symbols and functions have the following definitions: 


k 2 = 


h 2 = 


4 ar 


K(m ) = 
E{;m ) = 
Il(n, m) = 


£ 2 + (a + r) 2 ’ 
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>o yjl — m sin 2 6 

r /2 /-— 
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r /2 do 


(A3) 

(A4) 

(AS) 
(A6) 

(AT) 


(l — n sin 2 0) \/l — 77i sin 2 6 
The algorithms used by Mathematica to approximate elliptic integrals numerically use a 
completely different method than the adaptive integration algorithm of the NIntegrate 
function. Our comparison of Eqs. (J2J) and flSJ) with Eqs. (1A1|) and (1A2|) showed that both 


expressions gave relative errors of better than ICE 8 (eight significant digits) for points near 
and far from the solenoid’s surface. This accuracy was sufficient for the goals of this paper. 
Some further tests we made to validate our results were the following: 

1. We verihed that Eq. (jHJ) agrees accurately with the analytical expression for the mag¬ 
netic field on the axis of a continuous solenoid with end surfaces at z — ±L/2: 


B z {r = 0 ,z)= Hot x - 


z + L/2 


z-L/2 


(AS) 


y/(z + L/2) 2 + a 2 ^{z-L/2Y + a\ 

This expression includes the limiting case L —> 0 when the solenoid reduces to a current 
loop of radius a and again the numerical approximations to the continuous solenoid 
were found to be accurate in that limit. 

For z — 0, Eq. (jA8j) predicts that B z ~ Hot (1 — 2(a/L) 2 + ...) to lowest order 
in ( a/L ) 2 , i.e., the relative error (B Z (Q, 0) — Hot)/{Hot) decays as 1/L 2 for large L. 
This 1/L 2 convergence is similar to what is found for the solenoid’s external field. 


2. We verihed that Eqs. (121) and (|3ll ) converge to the corresponding components 


B r (r,z) = 


Hotn 


3 rz 


B z (r,z) = 


Hom 2 z 2 — r 2 


(A9) 


47T ( r 2 _ j _ - 2^ 5 / 2 ’ 2 ’ 47T ( r 2 _|_ - 2^) 5 / 2 

of the magnetic held of the continuous solenoid’s equivalent point magnetic dipole at 
its center r = z = 0, with magnetic moment 


m — 7ra 


! ) iL, 


(A10) 
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when the distance d = \/r 2 + z 2 of the point (r, z) of evaluation was large, d a, L, 
compared to the solenoid’s radius or length. That is, the solenoid’s external magnetic 
field correctly converges to that of a point magnetic dipole far from the solenoid, 
with a 1/d 3 decay in magnitude. Eq. (IA10|) was deduced by obtained by comparing 
Eq. (J6]) with the magnetic held of a point dipole po'm/ (47rr 3 ) for z = 0 in the far-held 
regime r L , a. 

3. We verihed that the total magnetic held B at points of certain symmetry had zero 
components, even though the numbers added to obtain the component were nonzero. 

We also verihed the convergence and accuracy of the 2D integral of magnetic energy, 
which presents two types of discontinuities, as shown in Fig. EH Firstly, the peak at r/a = 1 
and z = L /2 comes from a divergence of the radial component B r (r, z) of the magnetic held. 
Near the point of divergence, B r (r,±L/ 2) ss ln\r — a\, which does lead to a convergent 
value for the external magnetic energy. Since the divergence due to the zero-thickness of 
the continuous solenoid leads to a finite magnetic energy, the magnetic energy of a physical 
solenoid wound with wires of finite thickness will show results compatible to those presented 
here. 


Appendix B: Adding up radial components B r from different solenoids 

For both the ring geometry and hexagonal array, one has to add magnetic helds from 
different directions. We show the main issues for the ring geometry. The calculation of the 
total magnetic held B(p, 6 , z) at some point P = (p, 6 , z) inside the common cylindrical sur¬ 
face p = R of the ring of n solenoids described in Section HTTl provides a nice example of how 
to change information described in the cylindrical coordinate system (r, 0, z) of one system 
(a given solenoid) to information described by the cylindrical coordinate system (p, 9, z) as¬ 
sociated with the common cylindrical surface. The key issue is that the radial unit vector r* 
associated with the ith solenoid varies with i as shown in Fig. [THl so one has to accumulate 
the components of the vector B r i r j along the unit vectors p and 6 of the (p, 9, z) coordinate 
system. The ^-component B z is easier to compute since the z axes of the solenoids are all 
parallel to the z axis of the cylindrical surface. 

If the radius a n of each solenoid is given by Eq. El the axis of each solenoid lies on the 
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FIG. 16: Surface plot of the integrand (2'nr)[B 2 / of the magnetic energy, from 
Eq. © for the parameters a — 1, L — 4, and t — 1. The peak at r/a = 1 and z = L /2 
comes from an integrable divergence of the radial component B r (r,z ) of magnetic field (see 
Fig. [31(a)), while the line of discontinuity from z = 0 to z = L/2 along the line r = a comes 
from the ^-component B z , which changes sign discontinuously as one passes from the 
inside to the outside of the solenoid (see Fig. [31(b)). The two discontinuities in the energy 
integrand lead to a challenging integral to approximate numerically. 
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FIG. 17: Subdivided region of integration for calculation of external and internal energy. 
We use the symmetry of the solenoid to calculate the magnetic field only for the right side 

of the solenoid, z > 0. 
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FIG. 18: Illustration showing how the total magnetic held B tot at a point P = (p, 6, z ) 
inside the common cylindrical surface p = R with axis at point C is obtained by adding 
the vectors B r J ; = B r i (r t , z,)x % of the magnetic held created by the ith solenoid whose axis 
at point A has the coordinates (p, 6) = (R + a n , #*). The length r t is the radial coordinate 
of the point P with respect to the solenoid’s axis at A. 


circle p = R + a n and has angular coordinate 

2tt 

6i = —(i — 1), i = l,...,n. (Bl) 

n 

If we define the two-dimensional vectors 

p = (cos(0),sin(0)), (B2) 

6 — (— sin(0), cos(0)), (B3) 

P = PP , (B4) 

a i = (R + a n ) (cos (6 1 *), sin (^)), (B5) 


then the vector p —a* points from the axis A of the ith solenoid to the point P, its Euclidean 
length r t 

r i — \\p — a* ||, (B6) 

is the distance of the point P to the axis of the ith solenoid, and so the radial unit vector r t 
along the radial coordinate centered on the ith solenoid is given by 


With this notation, the radial and azimuthal components of the total magnetic held at the 
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point P are given by 


n 


B p — p ■ ^2 B r ,iXi, 

i =1 

(B8) 

n 

Be = 0'22 

(B9) 


i =1 


where B ri = B r (ri,z) is the radial component of the magnetic field Eq. (J2J) at a point a 
distance r = r* from the axis of the ith solenoid. 
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